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ABSTRACT 

We perform a linear magnetohydrodynamic perturbation analysis for a stratified magnetized envelope where 
the diffusion of heat is mediated by charged particles that are confined to flow along magnetic field lines. We 
identify an instability, the "coulomb bubble instability," which may be thought of as standard magnetosonic fast 
and slow waves, driven by the rapid diffusion of heat along the direction of the magnetic field. We calculate the 
growth rate and stability criteria for the coulomb bubble instability for various choices of equilibrium conditions. 
The coulomb bubble instability is intimately related to the photon bubble instability. The bulk thermodynamic 
| properties of both instability mechanisms are quite similar in that they require the timescale for heat to diffuse 

across a wavelength to be shorter than the corresponding wave-crossing time. Furthermore, over-stability occurs 
only as long as the driving resulting from the presence of the background heat flux can overcome diffusive Silk 
damping. However, the geometric and therefore mechanical properties of the coulomb bubble instability is the 
complete mirror opposite of the photon bubble instability. 

The coulomb bubble instability is most strongly driven for weakly magnetized atmospheres that are strongly 
convectively stable. We briefly discuss a possible application of astrophysical interest: diffusion of interstellar 
cosmic rays in the hot T ~ 10 6 K Galactic corona. We show that for commonly accepted values of the cosmic 
ray and gas pressure as well as its overall characteristic dimensions, the Galactic corona is in a marginal state 
of stability with respect to a cosmic ray coulomb bubble instability. The implication being that a cosmic ray 
coulomb bubble instability plays a role regulating both the pressure and transport properties of interstellar cosmic 
rays, while serving as a source of acoustic power above the galactic disk. 
Subject headings: MHD: instabilities - Galaxy: structure - ISM: cosmic rays 
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1. INTRODUCTION 



The idea that the diffusive flow of energy up through an atmosphere carries with it, the potential to de-stabilize acoustic motion, 
has been considered for quite some time (Baker & Kippenhahn 1962). The most familiar situation is where the flow of energy is 
transmitted by a radiative heat flux and the instability mechanism operates due to changes in the opacity that result from nearly 
O adiabatic changes in temperature and density, along the wave. The instability mechanism, known as the K-mechanism, is responsible 
,H , for the strong observed pulsations in RR-Lyrae and Cepheid variable stars. 

■ In the case where the instability is strong such that the pulsation amplitude is large, the overall structure and evolution of the 
*~ ' equilibrium flow may be significantly altered. If the saturation amplitude is small, information on the overall structure of the system 
J> . may be extracted. 

k> ' More recently, the inclusion of magnetic fields into perturbation analyses of optically thick, radiating, stratified flows has produced 
rN | some interesting results. In particular, accretion flows onto black holes and neutron stars are unstable to the photon bubble instability 
jj^ ■ (Arons 1992; Gammie 1998; Blaes & Socrates 2001). Furthermore, the photon bubble instability, which can be thought of as a 
standard magnetosonic wave that is driven over-stable by the presence of a stratified radiation field, may operate in equilibria that are 
weakly magnetized and/or highly sub-Eddington (Blaes & Socrates 2003). 

Balbus (2001) showed that incompressible Brunt- Vaisalla oscillations are unstable when horizontally-separated fluid elements 
are thermally-connected by coulomb conduction along field lines of relatively little dynamical importance. This "magneto-thermal 
instability" operates as long as the diffusion of heat along the field lines is rapid in comparison to the Brunt- Vaisalla frequency. For 
sufficiently weak fields, instability occurs for atmospheres with an outward decreasing temperature profile, rather than an outward 
decreasing entropy profile - as in the case of convective stability (see e.g. Parrish & Stone 2005; Chandran & Dennis 2006). 

In this work we explore the possibility that magnetoacoustic waves may be secularly driven when the envelope's flux of heat is 
mediated by anisotropic thermal diffusion that occurs solely along magnetic field lines. The stability of magnetized neutron stars, 
cosmic ray diffusion in the interstellar medium and the thermal structure of cluster-scale cooling flows are some examples of where 
the following instability analysis may be be relevant. 

The plan of this work is as follows. In ^2] we express our basic assumptions, write down the necessary conservation laws and 
then specify the parameters of the simple equilibria, which we then perturb. Furthermore, in §2.2l we place our analysis within the 
context of previous work on linear MHD theory where the effects of rapid heat diffusion was taken into account. In ^3] our linear 
analysis begins, with an emphasis on the thermodynamics of the coulomb bubble instability. In $4] our linear analysis continues with 
an emphasis on the mechanics of the coulomb bubble instability. In §4.41 we present a physical description of the coulomb bubble 
driving mechanism in juxtaposition with the photon bubble instability. In 5j3J we discuss the possibility that the Coulomb Bubble 
instability might play a role in determining the nature of interstellar cosmic ray diffusion. In 5j6]we summarize our results and discuss 
possibilities for further research. 
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2. ASSUMPTIONS, FUNDAMENTAL EQUATIONS AND PERTURBATIONS 

2.1. Conservation Laws and the Background 

We assume that the basic equations of ideal MHD apply. Also, the only source term in the first law of thermodynamics is provided 
by a diffusive heat flux. The conservation laws given below are the same as those found in Balbus (2001); they are 

^ + V-(pv) = (1) 

P^ + v-Vv^) =-VP+pg+i-(VxB)xB, (2) 

pW^+v-W) = -V Q, (3) 
Q = - X bb-V7\ (4) 

— = Vx(vxB), (5) 
at 

and 

V-B = 0. (6) 
Eqs. (H}-© enforce conservation of mass, momentum, and energy, respectively while eq. © enforces magnetic flux-freezing. 
Furthermore, eq. (|4|i constrains that the flow of heat follows a diffusion law along magnetic field lines and eq. I© enforces B to be 
purely solenoidal. Definitions of the various symbols in the expressions above, as well as others used throughout, are listed in Table 

m 

The form of the heat flux given by eq. is valid as long as two conditions are satisfied. First, the Larmor radius must be small 
in comparison to the mean free path of the particles endowed with relatively large amounts thermal mobility and therefore, Q || b. 
Second, the mean free path along B must be smaller than the characteristic scales of the problem e.g., the temperature scale height. 
Later on, when we consider perturbations, the mean free path along the B must be smaller than the wavelength of the perturbation in 
question. 

We assume that the background is in hydrostatic balance and the equilibrium magnetic field B is a constant. We have 

--VP = g, (7) 

P 

where changes in the gas pressure P are related to changes in density and temperature via an equation of state 

P=^L. (8) 

Throughout, we assume that the mean molecular weight \x is a constant so that changes in composition are ignored. We only consider 
a background that may be characterized as a "stellar envelope." That is, both the surface gravity and thermal energy flux throughout 
the medium is set to be a constant. In order to resemble a stellar envelope, the condition for radiative equilibrium must read 

V-Q = (9) 
and in addition, the mass of the atmosphere must be insignificant in comparison to the source of the gravitational field. 

Table 1 

Definitions of highly-used symbols 



Symbol Quantity 

p density 

v velocity 

P pressure 

T temperature 

s entropy per unit mass 

B magnetic flux density 

b magnetic field unit vector 

Q heat flux 
X thermal conductivity along field lines 

g gravitational acceleration 

uj wave frequency 

k wave vector 

Ci isothermal sound speed 

Va Alfven velocity 

k opacity 

. x(k-b) 2 

^.diff UJ + l l ,T(a s /dT) p 

OJ~ LJ 2 -(k-VA) 2 

4 Lagrangian displacement 

A = 8 + £ ■ V Lagrangian variation 

Aq asymptotic (in k) growthrate 
<5Q^ Eulerian heat flux perturbation _L B 
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2.2. Nature of Perturbations in Different Regimes 

Our background is simple and closely resembles a basic stellar atmosphere, threaded by an equilibrium magnetic field. The only 
major difference is that energy does not flow upwards via radiative diffusion, directly against the pull of gravity. Instead, the transfer 
of energy is mediated by charged particles of relatively high thermal mobility that drift along the equilibrium field. 

Now we consider dynamical fluctuations. To keep things simple, we examine local WKB perturbations whose time and spatial 
dependence are oc e I *' x ~" art ->, Furthermore, we only consider wavevectors that are two dimensional i.e., k = (k x ,k z ), where z lies in the 
vertical direction. 

2.2.1. Adiabatic Fluctuations: Alfven, Gravity and Magnetosonic Waves 

In the limit of slow thermal diffusion, linear perturbations to the first law of thermodynamics, eq. (01, take on a simple form 

-iujpT{5s+i-Vs) = -iujpTAs ~ 0. (10) 

The Eulerian component of As leads to the fluid's acoustic response for short-wavelength perturbations, while the Lagrangian 
component oc £ ■ V* leads to the incompressible gravity waves. 

Including magnetic forces, compressible (k ■ S\ 0) short-wavelength perturbations obey the magnetosonic dispersion relation 
given by 

~ 2 

to 2 = k 2 c 2 ^ + k 2 vl, (11) 

where c s is the adiabatic sound speed and all other symbols are defined in Table Q] The eigenvectors corresponding to the roots of 
the dispersion relation above are referred to as the fast and slow magnetosonic waves. 
Incompressible perturbations, with k • 5\ = 0, satisfy the dispersion relation 

k l 
k 2± 

where N\ v is the Brunt- Vaisalla frequency. Basically, the restoring force is provided by a combination of Alfvenic tension and 
buoyancy. 



^ 2 -7T^=0 (12) 



2.2.2. Highly Non-Adiabatic Incompressible Fluctuations: Balbus ' 2001 Analysis 

The fluctuations mentioned above are quite standard and do not therefore, merit further analysis or discussion. Interesting physical 
effects, such as damping and instability, occur when fluctuations with respect to the flow of matter and the flow of heat separate from 
one another i.e., when the perturbations become non-adiabatic. Somewhat counter-intuitively, non-adiabatic effects occur when the 
time scale for heat flow over a wavelength is shorter than the oscillation period of the perturbation in question. That is, the flow 
becomes non-adiabatic once the agent of heat transfer can either inject into, or remove energy from, the flow before the fluid has time 
to respond. 

Balbus (2001) studied incompressible waves (k ■ S\ = and c s — ► oo) in the limit where the diffusion time along the magnetic field 
is short in comparison to the wave crossing time. Another condition for Balbus' (2001) instability is that the Alfven time must be 
shorter than the wave crossing time so that Alfvenic tension cannot suppress the unstable growth of the perturbations. He found 
that Brunt- Vaisalla oscillations (or g-modes) are driven unstable if horizontally-separated fluid elements can quickly transfer heat 
from regions of relatively high to low temperatures, once they are perturbed in opposite directions along the background temperature 
gradient. The role of the magnetic field in this "magneto-thermal instability" is to serve as conduit of thermal energy between 
horizontally-separated fluid elements. 

2.2.3. Highly Non-Adiabatic Compressible Fluctuations: This work 

We now consider the manner in which compressible (slow and fast) magnetosonic waves are affected by the rapid flow of heat 
along magnetic field lines. In a sense, the instability and damping mechanisms covered in the next few sections can be viewed as 
a natural extension of either Balbus' (2001) or Blaes & Socrates' (2003) analysis of magnetoacoustic waves driven by the rapid 
diffusion of heat along the radiation pressure gradient. 

3. TOTAL PRESSURE PERTURBATION IN THE LIMIT OF RAPID HEAT CONDUCTION 

The rapid diffusion of heat alters standard magnetoacoustic motion through the pressure perturbation, which we write as 

»-(©*♦(£)-■ 

That is, we assume that the fluid's pressure is a function of of p and T only such that P = P(p,T). Furthermore, (dP/dpJ T = c 2 
where c, is the isothermal sound speed. For short- wavelength magnetosonic motion, Sp is constrained by the linearized expression 
for conservation of mass, while ST is constrained by the linearized first law of thermodynamics, which reads 

-iLunT(8s + £-Vs) = -ik-SQ, (14) 

where s and 5Q is the specific entropy and the perturbed coulombic heat flux, respectively. Later on, we show that the manner 
in which SQ responds to changes in magnetofluid variables ultimately determines the nature of the various driving and damping 
mechanisms. Below, we carefully discuss the form of SQ. 
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3.1. The Heat Flux Perturbation 
The equilibrium coulombic heat flux is given by 

Q = - X bb Vr=-X VT. (15) 

Here, X is the conductivity along field lines which lie in the direction R/B = b and X = X(x,b) is the thermal conductivity tensor. 
Upon perturbation, we have 



5Q = -SX-VT-X-V5T: 



— b (b • Vr) +5b (b • Vr) +b (5b ■ VT) + ib (k • b) ST 



(16) 



where the coulombic conductivity x is not taken to be a constant. We choose to parameterize the conductivity \ i n the following 

way, 

x = X() T"p m K P (17) 

where xo is a constant and n may be thought of as an "opacity" or in other words, a cross section per unit mass. Comparison with 
thermally diffusive damping and driving mechanisms for backgrounds in which radiation mediates the flow of energy is facilitated if 
we choose n = 3, m = -1, and p = -1 as in the case of radiative diffusion. In a way, we absorb the micro-physical differences (such 
as cross section) between radiative and coulombic diffusion into an opacity law - or a mean free path - for the particles that possess 
relatively large amounts of thermal mobility. 

From the above considerations, the expression for the conductivity perturbation is given by 



$P ^ K 
~f~~ 



(18) 



X P 1 K 

For now, we somewhat artificially set Sk = 0, which momentarily prevents us from studying a potential K-mechanism that results 
from coulomb diffusion. With this, the expression for the heat flux perturbation becomes 



SQ 



-jb(b- VT^j +5b(b- VT^j +b(sb- VT^j + ib (k ■ b) ST 



(19) 



Note that we have dropped the term in 5x oc 3 ST since it contributes at a lower order in (kH)~ l relative to the last term in the 
above expression, within the context of our WKB approximation. We do this with hindsight earned from Blaes & Socrates' (2003) 
analysis. The rapid diffusion (high-fc limit) of heat along temperature gradients necessarily implies that the temperature perturbation 
is relatively small by a factor of l/k in comparison perturbations of other magnetofluid quantities such as Sp, Sb and S\. 

3.2. The Temperature Perturbation 

By expanding eq. ([Pil l, the first law of thermodynamics, with the help of the expression for the heat flux perturbation SQ given by 
eq. dl9V we have 



-ujnT 



6p+ 



ds 
~dT 



8T + £-Vs 



= X 



p v 



k-b) (b- VT) + (k-6\i) (b ■ Vt) + (k ■ b) (sb ■ VT^j +i (k-b) ST{10) 
h quant 



It is useful to define a characteristic diffusion frequency, wc.diff, which quantifies the rate at which heat diffuses over a wavelength, 

2 

x(- 

<^C.diff = LO + l —rz T—s ■ (21) 



nT{ds/dT)- 



With this, we temperature perturbation becomes 



ST-- 



uj (ds/dp) T us £-Vs 



Wc.diff 



(ds/dT) 



-Sp 



X 



P 



k-b b-VT + (k-Sb ){b-VT )+ k-b ){Sb-VT 



^Cdiff 



{ds/dT) t 



disnT (ds/dT) 



(22) 



By knowing beforehand that the magnetosonic waves of interest can be treated to lowest order as standard fast and slow magnetoa- 
coustic waves tells us that u ~ v p ^k, where v p h is the phase velocity of the given wave. Therefore, in the rapidly diffusing limit, 



I^Cdiffl 



X k-b 



nT (ds/dT) 



> 0J. 



Now, the temperature perturbation reduces to 



ST 



npTuj ( ds\ Sp Sp 



X 



(23) 



(24) 



Note that the thermal coulomb conductivity \ was set to a constant in Balbus' (2001) analysis. 
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The form of the temperature perturbation given by eq. (f24l > is remarkably similar to the case of rapid radiative diffusion, which 
potentially leads to photon bubble-like phenomena i.e., 

npTuj ( ds\ Sp 5p 



57 ~T? 



(k-VT) 

X \OPj T P P 



RADIATIVE DIFFUSION. (25) 



In both cases, the first term which is oc (ds/dp) T leads to diffusive thermal damping. For the specialized case of radiative diffusion, 
this effect is known as Silk damping (Silk 1968; Weinberg 1971). The remaining terms in eq. d24l i and eq. d25T l - the one that contain 
gradients in temperature - lead to secular over-stable driving. For photon bubbles driven by rapid radiative diffusion, the term 
oc k- 'VTSp/p in the temperature perturbation represents "shadowing" or "pile up." More specifically, immediately downstream 
a density maximum, radiation piles up due to the local increase in extinction and since photon number is conserved, a deficit of 
radiation, or a shadow, occurs. Thus, the resulting radiation pressure differential across the density maxima allows for the possibility 
of work being performed on the fluid by the radiation field. In the case of anisotropic coulomb conductivity, a local change in the 
orientation of the magnetic field either permits or deters the flow of heat across a local density maximum. The orientation of Sb 
relative to constant density surfaces (that are _L k) ultimately determine whether or not the change in heat flow drives or damps the 
oscillation. 

It follows that the total pressure perturbation is divided into two parts; one which yields a standard acoustic response oc Sp while 
the other component oc ik~ l Sp leads to conductive driving and damping. That is, 

SP=SP :iC + SP (26) 

where 

SP. dc =(^f] Sp = cjSp and SP = (|^| ST. (27) 



4. COUPLING OF MAGNETOSONIC MOTION TO THE BACKGROUND COULOMB FLUX 

In what follows, we calculate the growth rates and stability criteria for the fast and slow coulomb bubble instability. We start 
with a brief overview of standard magnetosonic waves, where background gradients are ignored. We derive the basic properties and 
reveal the nature of the coulomb bubble instability by determining the ratio of the work done upon a magnetosonic wave by the 
driving mechanism to the wave energy (or wave action), which is equal to the ratio of the growth/damping rate to the oscillation 
frequency of the wave. We first arrive at this "work integral"-like ratio by examining the product of the Lagrangian pressure and 
density perturbation APAp, which closely resembles the -PdV work done on the wave. Then, we derive the work done on the wave 
through the quantity £ • <5f, where Sf is the linear driving force resulting from rapid coulomb diffusion. 

4.1. Basic Magnetosonic Waves 
Linearizing the continuity, momentum, and induction equations, as well as Gauss' Law, yields 

-iu;Sp+ip(k-S\) + S\- Vp = 0, (28) 

- iupS\ = -ik SP + g Sp + — (k x SB) x B , (29) 
4n 

-iLuSB = ikx(SvxB), (30) 

and 

k-(5B = 0, (31) 

respectively. 

Blaes and Socrates (2003) provide the mathematical apparatus for ascertaining the growth rate, stability criteria and basic physics 
of magnetosonic waves that are either driven or damped by the rapid diffusion of energy (see Socrates et al. 2005 and Turner et al. 
2005 as well). Solutions of the entire dispersion relation - for the case of radiative diffusion - show that in the high-A: limit, the growth 
rate Aq of the fast and slow wave, resulting from photon bubble driving, is a constant. Since the oscillation frequency for standard 
magnetosonic waves co oc k the ratio of the magnitude of the driving force to the magnitude of the magnetosonic restoring force is 
oc Ao/k oc 1 jk. It follows that the terms responsible for driving and damping in the linearized Euler equation are those that are oc ST, 
the buoyancy force oc gSp and the component of the Lagrangian density perturbation oc £ • Vlnp, since all of these quantities belong 
to linear forces that are l/k times smaller than the components responsible for the linear magnetosonic restoring force. In order to 
evaluate the growth rate, we insert the magnetosonic eigenvector into the O (k~ l ) linear forces, giving us the desired correction to the 
oscillation frequency. The technique outlined above is quite similar to the work integral approach in classic stellar pulsation theory 
(Unno 1989) for obtaining growth rates due to non-adiabatic driving of nearly adiabatic pulsations or to the method of calculating 
shifts in eigenfrequency in the linear perturbation theory of particles in quantum mechanics. 

We start by obtaining the form of the basic magnetosonic eigenvectors and their respective eigenfrequencies. In the high-£ limit, 
the dominant terms of the Euler equation are 



-icupS\ = — ik 



B SB 

5PK + 



k ■ B 

+ i SB, (32) 

4-7T 



4~ 

where SP- dc = cfSp measures the acoustic pressure response to a perturbation in density Sp. Furthermore, 

k-8\ B „ kB r 
Sp = p and SB = -(k-Sv) Sv, (33) 
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to leading order. With this, we write SB in terms of Sp and 5\ so that 



-ip — o\ ■ 



-ik 



I 2 , 2\ r 2^' v a) 



+ i'(k-v A )vA<5/0, 



(34) 



where va = B/y/4np is the Alfven velocity and Co 



~,2 



- (k • Va) 2 - In the short-wavelength limit, changes in temperature result- 



ing from coulomb conduction force the gas via a pressure gradient, which is || to k in the short wave-length limit. Therefore, 
the fluctuations of interest must possess some longitudinal, or compressible, component which automatically rules out purely in- 
compressible shear Alfven wave as a candidate for over-stable coulombic driving. We define a mode polarization e(u,k) such that 
S\ = e{uj , k)ip{u> , k) where ip is some complex-valued amplitude. We arbitrarily choose ip = Sp/p for the magnetosonic eigenvectors 



-cf + v A k-(k v A )v A 



where the two magnetosonic waves must satisfy the dispersion relation 

n 2 

2 ;2 2 W , ;2 2 

u> = k C: — - + k v A . 

In principle, the fast and slow waves posses components both || and _L to b = B/B = Va/va- 

4.2. The Magnetosonic Wave Equation Subject to Rapid Coulombic Diffusion: Asymptotic Growth Rate 



(35) 



(36) 



By taking the divergence of eq. (1291 , we immediately restrict our analysis to compressible short-wavelength perturbations i.e., the 
waves that can be driven unstable by the effects of thermal diffusion. We write 



;k • S\ = —ikr 



SP B SB 

p 4irp 



+ (k-g)^ 

P 



where we have made use of eq. (1311 1. We eliminate SB in favor of S\ with the help of eq. (13 Ot . which allows us to write 



-iui (k 5\) = -ik 



SP 2 (k-Sv) (k-v A ) 2 SP ,(g-v A )(k-v A )£p 
— + v, = ; = 



P 



+ (k-g)^. 

P 



(37) 



(38) 



By expanding the pressure perturbation into its acoustic SP &C and driving/damping SP contributions and by utilizing eq. fl28l ), we 
separate the standard magnetosonic wave equation from terms that may lead to driving and damping 



-e- 2 c 



; 2 2 

-lev. 



Sp 



~,2 



SP ,8y 



— ~ -iuoSy ■ Vlnp+/t < — Yiv A Vlnp-/ 



.(g-VA)(k-v A ) Sp 




(39) 



p uj or p 

The left hand side of eq. ( f39l > is the magnetosonic wave equation for a uniform background. The terms on the right hand side of 
eq. (|39l produce 0(k~ l ) corrections to the magnetosonic wave eigenvectors and oscillation frequencies as they are responsible for 
driving and damping. In order to evaluate the frequency correction arising from coulombic diffusion, we let u> — > ui +Aq, where Aq is 
the damping or driving rate that is independent of the wavenumber k. 

The condition for hydrostatic balance, eq. (0, allows us to recast the gravitational acceleration g in the following useful form 

I i / f)p\ 

g=-VP = cjV\np+- — VT. (40) 
P P\ dT Jp 

Together with the expression for the magnetosonic polarization vector given by eq. (l35l ). we eliminate all of the terms oc Vlnp in eq. 
d39l ). The perturbative prescription for the eigenvalues motivates conversion of the altered magnetosonic wave equation, given by eq. 
d39l ), into an expression for the damping or driving rate. After a bit of algebra, we have 



2A 



2uj z 



Sp 



uj 1 p 



.(k-v A ) 

-i — o — Va- 

UJ"p 



dP_ 

df 



VT 



+ i- 



dP_ 



VT 



S_p_ 
P ' 



By noting the form of the magnetosonic polarization vector given by eq. ( f35l >. we may write 



2A 



[2. 



^)] 



Sp 



k 2 cj 2 1 fdp 

dT 



uj 1 p 



AT 



(41) 



(42) 



where AT = ST + £- VT is the Lagrangian temperature perturbation. The expression for the growth rate given above may be written 
into a more a revealing form if we realize thatQ 

= ^ (43) 

which allows us to write 



Ao 



2 pS\-S\ 

4 Due to our choice of normalization (k ■ e) = ui, a relation that will be useful when we derive the growth rate later on. 



(44) 
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since 5\ = e5p/p and (k • = /Ap/p0 Eq. d44l i has a straightforward interpretation. The ratio of the growth rate to the oscillation 
frequency, Aq / u, is equal to the work done upon a nearly isothermal perturbation by changes in the flow's coulombic flux of energy, 
divided by the energy of that perturbation. 

Before we complete our calculation of the growth rate Aq, we recast the perturbation of the magnetic field unit vector in the 
following helpful forms 



<5b = bx 



SB A SB 
— x b = — 
BIB 



1-bb = 



<5B 



= ;(k-b)£- (l-bb) =- 



k 2 



(45) 



where eqs. d33l ) and d35l l were put to use. The Lagrangian temperature perturbation AT may be broken into a part that is responsible 
for diffusive Silk-like damping, which we denote as Ar damp and a portion A!T dnve which may lead to over-stable driving. We have 



Ar damp = Jr damp _ ■ n P TuJ f 5 _P. 



X 



d P 



(46) 



and 



Ar 



drive 



5T d ™ e + £- VT- 



k b 



-^(k-b)(b-VD + (k-(Sb)(b-vr) + (k-b) 



+£■ VT 



k b 



(k-0 k-b b-vr + k-b k- £• 1-bb b-Vr + k-b £■ 1-bb -VT 



+ £-VT 



6-vr] 



k b 



Ck-o (k-b)-(k-b) 2 (b^ 



Co 2 k 2 



b - vr )|.<5b 



Co 2 k 2 



£ 6Qj 



Sp/p 



AT d 



(47) 



X(k-b) Sp/p 

where we have made extensive use of eq. d45l l. In above expression for Ar dnve , the term <5Qj_ = SQ ■ ^1-bb^ is the component of 

the heat flux perturbation SQ that is perpendicular to the direction of the equilibrium field b (and thus Q). As eq. d47b implies, the 
perpendicular heat flux perturbation <5Q^ takes the form 



1-bb )-SQ = -x 1-bb 



5b lb ■ VT 



= -x(l-bbV i(k-b\£- (l-bb 



b- VT] =■ 



■x 



<5b(b-VrV (48) 



since by construction, 5b is _L to b. Thus, the form of Ar dnve tells us that work done upon a fluid element by SQ is proportional 
to the projection of the heat flux perpendicular to the equilibrium field, which is parallel to the perturbation of the magnetic field 
direction vector 5b, upon the fluid displacement Ultimately, the perpendicular heat flux <5Q^ must correspond to the linear driving 
force - a point that we further consider in the following section. 

Finally, we are in the position to write the growth rate Aq in terms of background quantities and the wave vector k of the oscillation 



iCj 2 k 2 



1 f dP 



pnTuj ( ds 



^ -2- 



(*■ 



VT 



2^[^-k 2 (cf + vl)}p\dT /p fcl 
and in terms of the background coulombic heat flux, the growth rate is given given by 



k b 



k b 



iu 2 k 2 



2uj[2uj 2 -k 2 (c 2 + v 2 )] \dT 



OP 



PX k-b 



pn Tuj 



as 

dp 



+ 2 — (k-Q) 



k 2 



k b 



(49) 



k 2 



(50) 



As previously mentioned, the first term, which contains the differential (ds/dp) T < 0, is responsible for diffusive Silk damping. The 
second term may contribute to either damping or driving, depending upon the direction of propagation relative to the vertical. Upon 
inspection, eq. (|49l indicates that for a given k, only one of the compressible magnetosonic waves is driven unstable, similar to the 
photon bubble instability. 

The relationship between the pressure and density perturbation of compressible MHD waves has allowed us to calculate the 
damping and growth rates arising from the action of rapid thermal energy transfer in a stratified background, resulting from coulomb 
diffusion along magnetic field lines. At the same time, we realized that the work done upon a fluid element by the background heat 
flux is proportional to the overlap between the fluid displacement £ and the component of the perturbed heat flux that is _L b. This 
clearly implicates <5Q^ as the driving force responsible for the over-stability. In what follows, we elaborate upon this point. 

5 Eq. 1441 applies for radiative diffusion as well. If ST is given by eq. {25}, then eq. )44t yields the photon bubble growth rate and stability criteria for the fast and 
slow magnetosonic waves. 



4.3. The Perpendicular Heat Flux, (5Qj_, as the Linear Driving Force and the Geometry of the Driving Mechanism 

Instead of examining the relationship between changes in pressure and volume for a driven magnetosonic wave, we consider 
the overlap of the driving force with the fluid motion itself. To isolate the linear driving force, it is convenient to work with the 
Lagrangian displacement The Euler equation reads 

0~-(k.v^) 2 ^-(cr + y2)k(k.e + VA(k-v A )(k.e + k(k.v A )(|-v A )-/k^^ ST -c 2 € • Vlnp^ - ig(k • © . (51) 

On the right hand side, the first four terms are responsible for the standard magnetosonic restoring forces while the other terms on 
the right hand side are responsible for secular driving and damping. Explicitly, 

(k- va) 2 £+ (c 2 + v 2 ) k (k • - va (k • v A ) (k © - k (k v A ) (© \ A ) = k 2 v 2 ph £, (52) 

where v 2 h = w 2 /k 2 is the magnetosonic phase velocity given by eq. (f36b . This allows us to write 



V + fc 2 v 2 h )£ 



dP 
P \df 



[kSr + (k-£)Vr] + /c 2 £x(kx Vln/j). 



(53) 



Note that the last term that is oc Vlnp is _L £ and therefore does not factor into the secular forcing of a given magnetosonic wave. It 
follows that 



(-c 2 + £ 2 v 2 h )^f 
where <5f is responsible for secular driving and damping and is given by 



Sf = - 



dP 
df 



[k57+(k-©VT]. 



(54) 



(55) 



As expected from eq. d44l ). forcing occurs as a result of changes in temperature. If we look at the projection of Sf along the fluid 
displacement we see that 



(56) 



P P \dTj p 

The connection between Sf and the Lagrangian pressure perturbation AP is clear and apparent. The left hand side of eq. d56*l l is 
the work done upon a fluid element irrespective of whether or not the motions involved are compressible or not. For the physics 
considered here, the driving force may be thought of as arising from changes in pressure due to changes in temperature, dictated by 
the flow of heat as prescribed by the first law of thermodynamics. In the previous section we found it useful to divide the Lagrangian 
pressure perturbation into a component solely responsible for diffusive Silk-like damping and one that may potentially lead to over- 
stable driving. Likewise, it is equally useful to decompose the secular driving force Sf into its analogous portions. With the help of 
eqs. d46i i and d4Tb we have 
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and 
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(58) 



The last two terms in the {} are equal to £ x (k x V7"), which is _L £ and cannot therefore, drive the fluctuations. Now, the linear 
driving force may be written as 
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(59) 



The form of <5f<jri V e in terms of the perturbed heat flux _L to b, given by <SQ = SQ ■ yl — bb J , allows us to plainly interpret the driving 

resulting from anisotropic coulomb diffusion. That is, <5fdrive results from the flux of energy that is _L to b, but projected along the 
wave vector k. 



Coulomb Bubbles 



9 



4.4. Essence of the Driving Mechanism and a Comparison with Photon Bubbles 

The growth rate for photon-bubble driving can identically be found by using eq. ( l44l i. but with an Eulerian temperature perturbation 
ST given by eq. ( f23T > rather than eq. ( l24l i. In a sense, the thermodynamics of compressible waves being driven by rapid thermal 
diffusion is similar in both cases. However, the mechanics and geometry of the respective driving mechanisms draws a definite 
distinction between the two instabilities. 

Previously, we compared the temperature perturbation in the limit of rapid coulomb diffusion to that of rapid radiative diffusion in 
order to distinguish between the physics of photon bubbles and of the coloumb bubble instability whose growth rate is given by eq. 
d49l ). For comparison, the photon bubble growth rate reads 




( +^(kxv A ).(kxF) 



(60) 



which is completely equivalent to the one-temperature MHD growth rates given by eqs. (93) and (107) of BS03. Note that the 
above expression is quite similar to the growthrate given by eq. ( [49) . However, the difference between the driving mechanism of the 
photon bubble instability and the coulomb bubble instability is clear upon examination of the driving component of the temperature 
perturbation for photon bubbles 

P SF 

Ar drive__i_i PHOTON BUBBLES. (61) 

X°PIP 

In this case, SF± = (1 - kk/A: 2 ) • <5F represents the component of radiative (heat) flux perturbation that is _L to the wave vector k. 

Figure s[T] and |2] display the geometry of the coulomb bubble instability and its relation to the photon bubble instability. Consider 
a fast magnetosonic wave in a stratified plasma where B 2 /8ir <C P in the limit of rapid anisotropic conduction. That is, the wave can 
be roughly thought of as a standard isothermal hydrodynamic sound wave with polarization vector e(ui,k) |j k with an oscillation 
frequency u ~ kc,. As this acoustic disturbance propagates throughout the atmosphere, the fluid motion along with the constraint of 
magnetic flux freezing induces a magnetic field perturbation SB. In general, SB possesses a component that is _L to the equilibrium 
field B, which then leads to a component of the heat flux SQ± that is _L b. It is this component of the linear heat flux SQ that 
is responsible for the over-stable driving mechanism of the coulomb bubble instability. As long as the component of the velocity 
perturbation S\ along k is endowed with a component that is _L to the equilibrium field, then driving may ensue in the event that Silk 
damping is overcome. 

Compare the geometry of the coulomb bubble mechanism described above with that of the photon bubble instability. Furthermore, 
consider Table 2, which compares the growth rates of both the coulomb bubble and photon bubble instability under various equilib- 
rium conditions. Apparently, the coulomb bubble instability may be thought of as the mirror image of the photon bubble instability. 
The prime mover for coulomb bubbles i.e., the physical quantity responsible for the driving, is the component of the heat flux pertur- 
bation 5Q± that is _L to the equilibrium field B. For photon bubbles, driving originates from the component of the radiative heat flux 
perturbation <5F^ that is _L to the wave vector k. Coulomb bubbles require that the component of the velocity perturbation S\ that is 
| to k possess some finite projection with the perpendicular heat flux SQ±, which is _L to B. For photon bubbles, the requirement 
is that a component of the velocity perturbation S\ that is || to 15 possess some finite projection with the perpendicular radiative heat 
flux <5F j_, which is _L to k. The most strongly-driven coulomb bubble wave is the fast wave in the limit where p ^> B 2 /87r (roughly a 
standard hydrodynamic sound wave) such that the velocity perturbation S\ is almost purely || to the wave vector k. Finally, the most 
strongly-driven photon bubble wave is the slow wave in the limit where p <C B 2 /87r such that the velocity perturbation S\ is almost 
purely || to equilibrium field B. 

5. AN EXAMPLE: INTERSTELLAR COSMIC RAY DIFFUSION 



Table 2 

Physical Properties and Conditions for Driving for both Coulomb Bubbles and Photon Bubbles." 

Mode Blanch Diffusion Law Plasma Beta Instability Criterion Asymptotic Growth Rate 
FAST Q = - X bb Vr p^>B 2 /&Tr Q>P c i (f-) 

SLOW 
FAST 



Q=- X bb-V7- p^B 2 /%* Q>[%\pvA 

Q = - X bb-V7- p<B 2 /&TT Q>PVA (t) (#) 

SLOW Q = - X bb-Vr p^B 2 /&7T Q>(X\pci (I) (f ) 

FAST F = - X Vr p>B 2 / 8w F ^\Jj^)p c i (j£ J (f ) 

SLOW F = - X Vr p>B-/8ir F> 



PVA 



FAST F = -x Vl p < B 2 /&tt F > p v A 

SLOW F=- X Vr p<B-/8ir F>pc; 



"These relations are approximate and only accurate to the order unity level. Terms that convey the geometric 
relationship between k, b, and the vertical gradient have been dropped. 
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fast wave 

B 2 /8ir < P 



FIG. 1 . — Geometry of the coulomb bubble instability. Here, the eigenfunction for a fast wave in a plasma with B 2 /8n/P < 1- nearly a standard hydrodynamic 
sound wave - is depicted. The velocity perturbation is perpendicular to surfaces of constant density and the perturbation of the perpendicular heat flux <5Q^ lies 
along 5b, the perturbation of field direction. Secular driving occurs because the wave in question possesses a finite projection along <5b, which is X to B. Thus, 
both the fast and slow waves are driven by the coulomb bubble mechanism whether or not the plasma is weakly magnetized or not. The most strongly driven case 
corresponds to the example above i.e., a fast wave in a weakly magnetized plasma. 
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FIG. 2. — Geometry of the photon bubble instability. Both the fast and slow waves are driven by photon bubble driving. The most strongly driven case - represented 
above - is the slow wave in a strongly magnetized plasma. 
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We now understand the inner workings of the coulomb bubble driving mechanism and the equilibrium conditions in which such 
driving can operate. Now, we put our analysis to use with the example of interstellar cosmic ray diffusion. 



5.1. Interstellar Cosmic Ray Diffusion and the "Two-Temperature" Approximation 

Throughout the analysis above, we assert that the thermodynamics of the coulomb bubble instability is identical to that of the 
photon bubble instability. Though this statement is valid, it is not entirely complete. Blaes & Socrates (2003) considered the possi- 
bility of photon bubble driving in two distinct thermodynamic regimes. If the micro-physical processes responsible for absorption 
and emission act quickly enough to maintain thermal equilibrium between the fluctuating radiation and gaseous particle field, then 
the thermodynamics of photon driving and damping are encapsulated by a common temperature. In other words, absorption and 
emission occur rapidly in comparison to the dynamical time so that the radiation and gas temperature are locked to one another. So 
far, our analysis of the coulomb bubble has taken place under this "one-temperature" approximation. 

If the particles endowed with a relatively high level of thermal mobility cannot come into close thermal contact with the bath 
of gaseous particles on a dynamical time, then the thermodynamics of radiative driving and damping is best described by a "two- 
temperature" approximation. An example of a system that is susceptible to photon bubble driving in the two-temperature limit is a 
magnetized envelope that is optically thick and whose only source of opacity is Thomson scattering. 

Interstellar cosmic ray diffusion is another example of where the transfer of a radiative energy is best described in a two- 
temperature approximation. That is, the Galactic distribution of cosmic rays cannot come into local thermal equilibrium with the 
distribution of gaseous particles due to the absence of absorption and emission processes between the two species. Therefore, the 
thermodynamic state of the combined matter + cosmic ray fluid is constrained by the first law of thermodynamics for the matter gas 
given by eq. ( |A3t in addition to eq. dA4b . the cosmic ray energy equation. 

Though cosmic rays cannot directly exchange energy with the matter gas, they certainly can exchange momentum via resonant 
scattering with magnetic irregularities on the Larmor scale. For ~ 1 — 10 GeV cosmic ray protons - the energy range that is responsible 
for the majority of the Galactic cosmic ray pressure - the Larmor radius r L is given by 



3 x 10 1 



(£ CR /GeV) 



(62) 



(B/pG) 

where £ CR and B is the cosmic ray energy and the large scale magnetic field strength, respectively. It follows that the length scale 
r L of the resonant magnetic irregularities are much smaller than the inferred energy-weighted cosmic ray mean free path A CR ~ 1 pc 
(Ginzburg et al. 1980; Strong & Moskalenko 1998) and the scale height of the Galactic corona, which is ~ a few kpc. 

The cosmic ray pressure perturbation 5P CR governs the exchange of momentum between fluctuations in the cosmic ray distribution 
and magnetofluid perturbations. The linearized Euler with the inclusion of the cosmic ray pressure force reads 



-iujpS\ = -ik [SP+5P CR ] + gSp+ — (k x SB) x B = -/k 
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(63) 



If Galactic cosmic rays were a "normal" radiation species that could come into close thermal contact the matter through absorption 
and emission processes, then local thermal equilibrium between the radiation and matter distribution may occur and then both 
radiation and matter distributions may be characterized by a common temperature T. Upon perturbation, both radiation and particle 
species share the same temperature perturbation ST as long as emission and absorption processes act quickly on dynamical timescales 
of interest (Blaes & Socrates 2003). For interstellar cosmic ray diffusion, the absence of any absorption and emission opacity 
necessarily implies that cosmic ray pressure perturbation SP CR and gas temperature perturbation ST are separate and decoupled. 

Below, we discuss the possibility that the diffusion of interstellar cosmic rays in the Galactic corona leads to a "two-temperature" 
version of the coulomb bubble instability. 



5.2. Pressure Perturbation, Stability Criteria and Growth Rate 

The diffusion of interstellar cosmic rays in the galactic disk and halo may lead to the two-temperature analogue of the photon 
bubble instability. In this case, cosmic rays, rather than photons, are responsible for the diffusive radiative transfer of energy. The 
action of cosmic rays scattering off of resonant magnetic irregularities leads to an "opacity" responsible for mediating momentum 
exchange between the radiation species and the gas. In the photon case, the relevant interaction is governed by Thomson scattering 
of photons with electrons. Rather than the radiation field diffusing along the gradient in the radiation pressure, cosmic rays diffuse 
along the projection of the cosmic ray pressure gradient that coincides with the equilibrium magnetic field. In what follows, we 
produce an abbreviated derivation of the growth rate Aq for the coulomb bubble instability resulting from cosmic ray diffusion in the 
galactic corona. 

The perturbation of the cosmic ray conductivity is given by 

Sx CR = Sp d\nn CR Sp 
X C R P 9\np p ' 

In Appendix [Al we take note that <91nK CR /91np ~ — 1 for the Galactic halo. The implication being that Sx CR — 0. However, we 
maintain the parameterization of Sx CR given by eq. d64l i in order to separate the effects of a "cosmic ray K-mechanism" and coulomb 
bubbles^ 

6 Note that the physical dimensions of the thermal conductivity \ and Xcr are different. Technically, x C r i s a diffusivity rather than a conductivity. 
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The perturbation to the cosmic ray energy flux may be written as 



^QcR = ~Xc 
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(65) 



which then allows us to determine the cosmic ray pressure perturbation <5P CR in terms of Sp and <5b. The linear cosmic ray energy 
equation reads 



- iu5P a + A -iuP c , j-'uii- (VP CR - 1 ^ Vp 



-zk -<5Q C 



(66) 



The second term on the left hand side in above expression leads to diffusive Silk damping, while the third term is responsible for 
Brunt- Vaisalla oscillations. In jj3.2| we took the high-/c limit of rapid diffusion of the first law of thermodynamics in order to isolate 
the temperature perturbation 5T in terms of the density and magnetic field unit vector perturbation 5p and 5b, respectively. In order to 
study coulomb bubbles in the "two-temperature" approximation appropriate for cosmic ray diffusion, we likewise isolate the cosmic 
ray pressure perturbation <5P CR in terms of Sp and 5b in the limit of rapid diffusion. By combining eqs. d65l l and d66*l l. the cosmic ray 
pressure perturbation becomes 
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where o;cR,diff = ^ + 'Xcr ( ^ ' ^ ) * s tne characteristic cosmic ray diffusion frequency. In the limit of rapid diffusion, we take 



(67) 
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Above, we take advantage of the fact that to leading order, the frequency of an acoustic oscillation is oc k. With this, we have 
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(68) 



(69) 



The similarity is evident between <5P CR above and the two-temperature radiation pressure perturbation 6P m a given by Blaes & Socrates 
(2003) 



5P^\8E = -i- 2 
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3 »P P 

where kf is the flux mean opacity as defined in Blaes & Socrates (2003). Note that the Silk damping term of eq. d69l is smaller than 
the one found in eq. d70T i by a factor of three, which results from our choice of parameterization of the cosmic ray diffusivity x CR - 
We are now in the position to calculate the growth rate and stability criteria for the cosmic ray coulomb bubble instability. Following 
the analysis that led to eq. (l44l i. the ratio of the asymptotic growth rate to acoustic oscillation frequency becomes 

1 — A ^CR 

-. (71) 

2 p5\ ■ 5v 

Note that in the two-temperature case the velocity eigenvector oscillates at magnetoacoustic frequencies that are the solution to eq. 
( fTTT i. The quantitative difference is that the isothermal gas sound speed c, is replaced with the adiabatic gas sound speed c s , since the 
gas neither generates nor loses any heat itself, the gas temperature perturbation then contributes to the acoustic response. In terms of 
its individual constituents, the growth rate may be written as 
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(72) 



The form of A ( > not only serves as an asymptotic growth rate, but as a stability criteria as well. The first term in the {} leads to 
Silk damping, the second term results in the coulomb bubble instability and the last term is responsible for a K CR -mechanism. The 

quantity 5Q CR .± = ^1 - bb^ • <5Q CR is the cosmic ray flux perturbation that is _L to the equilibrium field B (and equilibrium cosmic ray 

flux Q CR ). It follows that the geometry and mechanics of the two-temperature cosmic ray coulomb bubble instability is identical to 
the one-temperature coulomb bubble instability. 

In terms of the wave vector and specified equilibrium parameters, the asymptotic growth rate becomes 
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5.3. Marginal Stability of the Galactic Corona 

As stated in Appendix[A] Galactic cosmic rays almost uniformly fill a halo (or corona) with a characteristic thickness ~ a few kpc. 
Over these scales, the coronal gas is hot (T ~ 10 6 K) and tenuous « ~ 10~ 3 cm~ 3 = 773 and the corresponding adiabatic gas sound 

speed c s ~ 10 7 r 6 '^ 2 cm/s, where T(, is temperature in units of 10 6 K. For our stability analysis, we take the halo equilibrium field to 

be uniform with a characteristic value of B ~ 1/iG = B^c, corresponding to an Alfven velocity va < 10 7 B AI G«3 1,/2 cm/s, somewhat 
less than the adiabatic sound speed c s . The average local pressure for gas and cosmic ray protons in the interstellar medium is 
roughly ~ leV/cm 3 ~ 10~ 12 erg/cm 3 . Due to the large ~ a few kpc scale height for the cosmic rays, this value of pressure also 
corresponds to the coronal value of the cosmic ray pressure P CR ~ leV/cm 3 . Compare this with the value of the coronal gas pressure 
P ~ 10 _1 «3 TeeV/cm 3 . Thus, the ratio of cosmic ray to gas pressure is P CR /P ~ 10 in the kpc-scale galactic corona, a somewhat 
surprising value. In short, the Galactic corona is a cosmic ray pressure supported atmosphere. 

In the diffusion approximation, the cosmic ray flux Q CR for a plane-parallel geometry satisfies the relation 

Z f L Qc~^- (74) 

For a disk-like geometry the corresponding cosmic ray luminosity L CR is roughly 

L CR ~5x lO^PjT^oergs- 1 (75) 

where T3.5, Pi and^io is the cosmic ray optical depth in units of 10 3 5 , halo cosmic ray pressure in units of 1 eV/cm 3 and characteristic 
cylindrical emitting radius in units of 10 kpc, respectively. We now have the necessary ingredients to determine whether or not the 
kpc-scale galactic corona is over-stable to either the coulomb bubble instability or the K-mechanism. When calculating the growth or 
damping rates resulting from cosmic ray diffusion, we assume that the equilibrium field is relatively uniform over a wavelength and 
that it does not contribute to hydrostatic balance as well. 

In j ^4.4| we made note of the fact that the coulomb bubble instability is most strongly driven for waves that closely resemble simple 
propagating hydrodynamic sound waves i.e., the fast wave in the limit where c 2 s 3> v\. In this case, the instability condition for 
cosmic ray coulomb bubble driving from eq. (l73l is approximately 

Q CR >P CR c s -^ — >c s . (76) 

T 
1 CR 

In the above expression we ignore geometrical and other factors of order unity. Interestingly, the quantity c/r CR , which represents the 
diffusive drift velocity vd over the coronal scale height, is close in value to the adiabatic gas sound speed i.e., vq ~ 10 7 cm/s ~ 

c s r~j 10 7 r 6 1,/2 cm/s. Altogether the condition for instability is satisfied as long as T^ 2 t$.s < 1. The growth rate r FAST for the fast wave 
coulomb bubble instability in the c 2 s ^S> v\ limit is given by 

If ~ 3 x 10-7 V^ir VmH ^ yrs_1 ' (77) 

where V300 is the circular velocity in units of 300km/s and Ht, is the coronal scale height in units of 3 kpc. The growth rate r FAST 
is dynamical and indicates that the fast wave coulomb bubble instability can transform a low amplitude propagating fast wave, with 
phase velocity equal to c s , into a large amplitude sonic disturbance over a relatively few wave crossing times. 

5.4. Coulomb Bubbles vs. then-Mechanism 

The approximate stability criteria and growth rate for a k cr - mechanism is identical to that of the cosmic ray coulomb bubble 
instability in the corona of the Galaxy. Fast waves in the weakly magnetized limit, with u ~ kc s , are the strongest growing excitations 
for both cases. However, eq. d73b informs us that upward propagating fast waves are driven by the coulomb bubble instability while 
downward propagating fast waves are driven by the K CR -mechanism, for <91nK CR / <91np ~ — 1, 

At the base of the corona, the Galactic disk is highly inhomogeneous and turbulent. The turbulence in the interstellar medium pri- 
marily results from explosive phenomena - expanding HII regions, line-driven stellar winds and core-collapse SNe - that originates, 
in some way, from massive stars. There are ~ 10 4 O stars in the Milky Way within a galactocentric radius of ~ 10 kpc. Therefore, 
the characteristic separation between luminous massive stars ~ 100 pc is roughly the scale lj at which the interstellar turbulence is 
stirredQ 

With respect to the Galactic corona, the dense multiphase disk serves as a source of acoustic radiation that can then be amplified 
in the corona by either the cosmic ray coulomb bubble instability or the K CR -mechanism. The direction of propagation for the seed 
k~ l ~ I T ~ 100 pc-scale acoustic fluctuations is primarily upwards with respect to gravity in the Galactic corona since the dense disk 
lies at its base. It follows that the cosmic ray coulomb bubbles are the the relevant mechanism that amplifies propagating acoustic 
perturbations arising from the Galactic disk rather than a K CR -mechanism. 

6. SUMMARY 

We have identified a new acoustic over-stability, which we refer to as the "coulomb bubble instability," that is driven by the rapid 
diffusion of energy along magnetic field lines. From a linear instability analysis, we calculate the condition for over-stability and 
the growth rate of the coulomb bubble instability. Driving occurs when the ratio of heat flux to entropy is relatively large for the 
equilibrium, implying that the heat flux is the physical quantity that provides the free energy required for the over-stability. The 

7 This approximation for It is valid as long as the lifetime of massive stars is short in comparison to the dynamical timescales of interest. 
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coulomb bubble instability may be thought of as a standard magnetoacoustic wave that is driven by SQ± = ^1-bbJ ■ SQ, the linear 

perturbation of the coulomb heat flux that is _L to the direction of the equilibrium field b = B/B (and thus Q = -xbb ■ VT as well). 
The growth rate of the coulomb instability is dynamical for the fast wave in the B 2 /8ir <C P limit - roughly a standard hydrodynamic 
sound wave - such that the growth time is roughly the sound crossing time over a gas pressure scale height. 

The properties of the coulomb bubble instability are, at the same time, strikingly similar and starkly different from the photon 
bubble instability. The thermodynamic properties of both instabilities are identical in that the rapid diffusion of thermal energy in 
a stratified flow is needed in order for driving to occur. The work done upon the fluid by the driving force for both instabilities 
is proportional to the product of the Lagrangian density and temperature perturbation ApA7\ Furthermore, and like the photon 
bubble instability, both the fast and the slow wave are susceptible to over-stability from the coulomb bubble mechanism irrespective 
of whether or not the background flow is strongly magnetized. 

Work is done upon the wave by the coulomb bubble linear driving force provided by <5Qj_, the component of the heat flux 
perturbation that is _L b. The fast and slow magnetosonic waves are polarized along both the wavevector k due to compression, and 
along the direction of the magnetic field b, which results from magnetic tension. In general, SQ± has some finite projection that 
is coincident with the wavevector k - as long as k |f b - which allows driving to occur, since all compressible fluctuations exhibit 
motion that is in part || k. On the other hand, the photon bubble linear driving force originates from <5F^ = (l-kk/£ 2 ) ■ <5F, the 
component of the radiative heat flux perturbation that is _L k. Since the magnetic tension force partially polarizes magnetosonic 
waves along b, the radiative driving force due to 6¥± is coincident with the motion of the fluid, which then allows the radiation field 
to perform work upon the fluid. Clearly, in describing the geometry of the driving mechanisms for the coulomb bubble and photon 
bubble instability, the vectors k and b are interchanged with one another in every possible way. From this, we conclude that the 
coulomb bubble instability is the mirror opposite of the photon bubble instability. 

Both the fast and slow magnetosonic waves are driven over-stable by the coulomb bubble mechanism. The fast mode in the weakly 
magnetized limit is the most strongly driven case, with a growth rate ~ g/c g , where c g is the gas sound speedQ If the equilibria in 
question resembles a stellar envelope, then both the coulomb bubble instability and Balbus' magnetothermal instability possess 
comparable growth rates. However, if the particles responsible for the transfer of thermal energy along field lines (e.g. cosmic rays) 
are not the particles responsible for providing the acoustic response (the gas), then the growth rate of the coulomb bubble instability is 
larger than the magnetothermal instability by a factor that is <~ c CR /c, ~ \/P CR JP, for the case of cosmic ray diffusion in the Galactic 
halo. 

The coulomb bubble instability may thrive in a wide variety of astrophysical environments. For example, it is possible that coulomb 
bubbles play an important role in the modifying the thermal transport properties of the envelopes of magnetized neutron stars, central 
cooling flows of galaxy clusters and interstellar cosmic ray diffusion. In the latter case, we show that the kpc-scale galactic corona is 
a setting in which a cosmic ray coulomb bubble instability possibly operates. For commonly accepted values of cosmic ray and gas 
pressure, cosmic ray luminosity and coronal gas temperature, the Galactic corona is in a state of marginal stability with respect to a 
cosmic ray coulomb bubble instability. The implication being that cosmic ray coulomb bubbles self -regulates the interstellar cosmic 
ray pressure on Galactic scales. 

6.1. What Is To Be Done? 

Our local linear analysis is only the first step in realizing whether or not the coulomb bubble instability has practical applications. 
At this simple exploratory level, we only examine local propagating WKB waves at linear order. The next step is to perform a global 
linear analysis, where a linear mode is really the combination of an upward and downward propagating wave, whose properties are 
further augmented by the specified boundary conditions. Of course, linear theory can only tell us whether an instability exists. A 
MHD computer algorithm that incorporates the effects of rapid anisotropic diffusion in a stratified atmosphere is most likely the 
correct tool in fully studying the non-linear outcome of the coulomb bubble instability. 
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APPENDIX 

AN ENVELOPE PARTIALLY SUPPORTED BY COSMIC RAYS 

In order to describe the physics of cosmic ray diffusion in a stratified atmosphere, we amend the conservation laws and equilibrium 
conditions of $2] Conservation of momentum now reads 

d\ \ 1 

_+vVvj =-V(P+P CR )+ P g+ — (VxB)xB, (Al) 

where P CR is the cosmic ray pressure - dominated by ~ 1 - 10 GeV cosmic ray protons in the interstellar medium. The first law of 
thermodynamics for the gas 

pT S +V ' VS ) = ° (A2) 
now implies that fluid-matter perturbations are adiabatic. Evolution of the cosmic ray pressure P CR is further constrained by a cosmic 
ray-energy equation 

^+vVP cr + 1p cr V-v = -V-Q cr , (A3) 
where we assume that adiabatic index of the cosmic rays is that of a relativistic ideal gas i.e., 7 CR = 4/3. The cosmic ray energy flux 

Q C r is g iven b y 

Q CR =- XcR bb-VP CR , (A4) 

which is similar in form to the coulombic heat flux given by eq. (0|. 

The cosmic ray conductivity x CR deserves discussion as its prescription for the Galaxy differs from the prescription for thermal 
conductivity \ described in 33.11 Measurements of the cosmic ray transport in the Milky Way indicate a cosmic ray mean free path 
A CR ~ lpc. Interestingly, models of Galactic cosmic ray transport are consistent with A CR roughly being a constant over a "halo" scale 
height of ~ a few kpc. It follows that the cosmic ray "opacity" k cr must vary inversely with density i.e., k cr oc p . The immediate 
suggestion is that a "cosmic ray K-mechanism" may operate in the galactic halo since d\nn CR / d\np ^0 and is given by 

d\nn ro 

-=^^-1. (A5) 
amp 

Finally, the only other basic assumption that differs from those given in ^involves hydrostatic balance 

--V(P+P CR )=g. (A6) 
p 

That is, partial support of the atmosphere is given by an equilibrium cosmic ray pressure gradient. 



